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Abstract 

Recent years have seen extensive applications of the Nambu-Jona-Lasinio (NJL) model in the 
study of matter at high density. There is a good deal of interest in the predictions of diquark 
condensation and color superconductivity, with suggested applications to the study the properties 
of neutron stars. As the researchers in this field note, the NJL model does not describe confinement, 
so that one is limited to the study of the deconfined phase, which may set in at several times nuclear 
matter density. Recently, we have extended the NJL model to include a covariant confinement 
model. Our model may be used to study the properties of the full range of light mesons, including 
their radial excitations, in the 1-3 GeV energy domain. Most recently we have used our extended 
model to provide an excellent fit to the properties of the 77(547) and r/(958) mesons and their radial 
excitations. The mixing angles and decay constants are given successfully in our model. In the 
present work our goal is to include a phenomenological model of deconfinement at finite matter 
density, using some analogy to what is known concerning "string breaking" and deconfinement at 
finite temperature. Various models may be used, but for this work we choose a specific model for the 
density dependence of the parameters of our confining interaction. We perform relativistic random- 
phase-approximation (RPA) calculations of the properties of the 7r(138), K(495), /o(980), ao(980) 
and i^g(1430) mesons and their radial excitations. In the model chosen for this work, there are no 
mesonic states beyond about 2p^M, where pnm is the density of nuclear matter. (The density for 
deconfinement in our model may be moved to higher values by the change of one of the parameters 
of the model.) This inability of the model to support hadronic excitations at large values of the 
density is taken as a signal of deconfinement. In addition to the density dependence of the confining 
interaction, we use the density-dependent quark mass values obtained in either the SU(2) or SU(3)- 
flavor versions of the NJL model. We stress that other assumptions for the density dependence 
of the confinement potential, other than that used in this work, maybe considered in future work, 
particularly if we are able to obtain further insight in the dynamics of deconfinement at finite 
matter density. 

PACS numbers: 12.39.Fe, 12.38.Aw, 14.65.Bt 
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I. INTRODUCTION 



In recent years we have developed a generalized Nambu-Jona-Lasinio (NJL) model that 
incorporates a covariant model of confinement [1-5] . The Lagrangian of the model is 

C = q(i$> - m°)q + ^ ^[(gA'g) 2 + {qi^Xqf} 

i=0 

-^E[(^7^) 2 + (^757,g) 2 ] 

i=0 

+ ^{det[g(l + 75 )g] + det[g(l - 75)?]} 

+£conf , (1-1) 

where the = 0,---,8) are the Gell-Mann matrices, with A = a/2/3 1, m° = 

diag (mj, m° d1 m° s ) is a matrix of current quark masses and C con f denotes our model of con- 
finement. Many applications have been made in the study of light meson spectra, decay 
constants, and mixing angles. In the present work we extend our model to include a de- 
scription of deconfinement at finite density. 

There has been extensive application of the NJL model in the study of matter at high 
density, with particular interest in diquark condensation and color superconductivity [6-9]. 
These studies find application in the study of neutron stars. The NJL model is the model of 
choice, since little insight into the properties of matter at finite density can be obtained in 
lattice simulations of QCD. This problem is associated with the introduction of a chemical 
potential, which makes the Euclidean-space fermion determinant complex. 

The use of the NJL model in the hadronic phase of matter is limited, since the standard 
version of the model does not contain a model of confinement [10-12]. It is clearly of value 
to extend the NJL model so that one can study the full range of densities of interest at 
this point in time. We are encouraged in this program by recent results, obtained in lattice 
simulations of QCD with dynamical quarks, that provide information on the temperature 
dependence of the confining interaction [13]. It is generally believed that the presence of 
matter will play a role similar to that of finite temperature, with deconfinement taking place 
at some finite density, which might be several times that of nuclear matter. In the present 
work we make a specific assumption concerning the density dependence of the confining 
field and then calculate meson spectra in the presence of our density-dependent confining 
interaction. We also take into account the density dependence of the constituent quark 
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masses, which is calculated in the SU(2) or SU(3)-flavor version of the NJL model. As 
is well known, the presence of matter leads to a reduction in the magnitude of the quark 
vacuum condensates, which represents a partial restoration of chiral symmetry in matter. 

Our calculations of the properties of mesons in matter is made using a covariant random- 
phase-approximation (RPA) formalism, which we have developed for the study of mesons in 
vacuum [1-5]. The organization of our work is as follows. In Section II we provide a short 
review of our treatment of Lorentz-vector confinement in our generalized NJL model. In 
Section III we describe the variation of the up, down and strange quark constituent quark 
masses in matter. In Section IV we discuss some recent work concerning the temperature 
dependence of the confining interaction, as obtained in lattice simulations of QCD with dy- 
namical quarks. We also specify the density dependence of the confining field that we use in 
this work in Section IV. In Section V we comment upon the phenomenon of pion condensa- 
tion. (In our work we introduce a small density dependence of the coupling constants of the 
NJL model to simulate effects that prevent the formation of a pion condensate in nuclear 
matter.) In Section VI we discuss our covariant RPA calculations of meson properties in 
vacuum and indicate how these calculations are modified in matter. Results of our RPA 
calculations of the properties of pseudoscalar mesons in matter are presented in Section VII, 
while Section VIII contains similar results for scalar mesons. In the case of scalar mesons, 
we study the a (980), / (980), and Kq(1A30) mesons and their radial excitations. Finally, 
Section IX contains some further discussion and conclusions. 



II. MODELS OF CONFINEMENT 



There are several models of confinement in use. One approach is particularly suited to 
Euclidean-space calculations of hadron properties. In that case one constructs a model of 
the quark propagator by solving the Schwinger-Dyson equation. By appropriate choice of 
the interaction one can construct a propagator that has no on-mass-shell poles when the 
propagator is continued into Minkowski space. Such calculations have recently been reviewed 
by Roberts and Schmidt [14]. In the past, we have performed calculations of the quark and 
gluon propagators in Euclidean space and in Minkowski space. These calculations give rise 
to propagators which did not have on-mass-shell poles [15-18]. However, for our studies of 
meson spectra, which included descriptions of radial excitations, we found it useful to work 
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in Minkowski space. 

The construction of our covariant confinement model has been described in a number of 
works [1-5]. We have made use of Lorentz- vector confinement, so that the Lagrangian of our 
model exhibits chiral symmetry. We begin with the form V c (r) = nrexp[— fir] and obtain 
the momentum-space potential via Fourier transformation. Thus, 



V c (k-k') 
with the matrix form 



-8ttk 



[(k- k') 2 + fi 2 } 2 {(k-k') 2 + fi 



213 



(2.1) 



V c (k-k')=Y(W C (k-k'h»(2), 



(2.2) 



appropriate to Lorentz- vector confinement. The potential of Eq. (2.1) is used in the meson 
rest frame. We may write a covariant version of V c (k — k') by introducing the four- vectors 

(k ■ P)P fl 



k» = k»- 



P 2 



(2.3) 



and 



(k' ■ P)P» 
P 2 



(2.4) 



Thus, we have 

V c (k-k') = -8tik 



4/i 2 



(2.5) 



[—(k - k') 2 + /i 2 ] 2 [—(k - k 1 ) 2 + jj, 2 ] 3 

Originally, the parameter /x = 0.010 GeV was introduced to simplify our momentum-space 
calculations. However, in the light of the following discussion, we can remark that \i may 
be interpreted as describing screening effects as they affect the confining potential [13]. 

In our work, we found that the use of k — 0.055 GeV 2 gave very good results for meson 
spectra. Here, k for the Lorentz-vector potential is about one-fourth of the value of k for 
Lorentz-scalar confinement. This difference arises since the Dirac matrices 7^(1)7^(2) in 
Eq. (2.2) give rise to a factor of 4 upon forming various Dirac trace operations, so that the 
effective value of the string tension is about the same in both Lorentz-scalar and Lorentz- 
vector models of confinement. 

The potential V c {r) = /trexp[— fir] has a maximum at r = 1/fj,, at which point the value 
is Vmax — n/pte — 2.023 GeV. If we consider pseudoscalar mesons, which have L = 0, the 
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FIG. 1: a) Bound states in the confining field (wavy line) may be found by solving the equation 
for the vertex shown in this figure, b) Effects of both the confining field and the short-range NJL 
interaction (filled circle) are included when solving for the vertex shown in this figure. 

continuum of the model starts at E cont = mi + m 2 + V max , so that for mi = m 2 = m u = 
m d = 0.364 GeV, E cont = 2.751 GeV. It is also worth noting that the potential goes to zero 
for very large r. Thus, there are scattering states whose lowest energy would be mi + m 2 . 
However, barrier penetration plays no role in our work. The bound states in the interior of 
the potential do not communicate with these scattering states to any significant degree. It 
is not difficult to construct a computer program that picks out the bound states from all 
the states found upon diagonalizing the random-phase-approximation Hamiltonian. 

Bound states in the confining field may be found by solving the equation for the mesonic 
vertex function shown in Fig. la. Inclusion of the short-range NJL interaction leads to an 
equation for the vertex shown in Fig. lb. We will return to a consideration of Fig. lb when 
we discuss our covariant RPA formalism in Section VI. 

III. CALCULATION OF CONSTITUENT QUARK MASS VALUES 

In this Section we report upon our calculation of the density dependence of the constituent 
quark masses of the up (or down) and strange quarks. The role of confinement in the 
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calculation of the constituent mass was studied in an earlier work in which calculations were 
made in Euclidean space [19]. The results were similar to those obtained in Minkowski-space 
calculations in which confinement was neglected and it is the latter calculations which we 
discuss here. 

The equations for the quark masses in the SU(3)-flavor NJL model are [11] 



m u = m" — 2Gs(uu) — Go(dd){ss) , 
m d = m° — 2Gs(dd) — Gd(uu)(ss) , 
m s = m° — 2Gs(ss) — Gd(uu) (dd) , 



(3.1) 
(3.2) 
(3.3) 



where (uu), (dd) and (ss) are the quark vacuum condensates. For example, with N c = 3, 

d 4 k 



(3.4) 



r) 4 k 2 — m\ + ie 

If this integral is evaluated in a Minkowski-space calculation, a cutoff is used such that 
\k\ < A 3 . Thus, 



, ,r d3k 



) 3 2E u (k) ' 



(3.5) 



etc. Here E u (k) 



k 2 + m 2 



1/2 



For studies at finite density, we consider the presence of two Fermi seas of up and down 
quarks with Fermi momentum kp. We also take m° u = m° d and obtain the density-dependent 
equations, with (uu) p = (dd) p , 



m u (p) =m° u - 2G s (uu) p - G D (dd) p (ss) p , 
m s (p) =m° s - 2G s (ss) p - G D (uu) p (dd) p . 



(3.6) 
(3.7) 



Equation (3.5) is now replaced by 



(uu) f 



-AN r 



1/2 



As d 3 k m u (p) f kF d 3 k m u (p) 



o (2vr) 3 2E u (k) Jo (2vr) 3 2E u {k) 



(3.8) 



with E u {k) = k 2 + m„(p) . On the other hand, since we do not consider a background 
of strange matter, we have 

" As d 3 k m s {p) 



{ss) p = -4iV c 



r 3 _^ 

Jo (2tt 



) 3 2E s (k) ' 



(3.9) 



with E s (k) = k 2 + m 2 s (p) 

We may argue that, with respect to our mean-field analysis, the Fermi seas of up and 
down quarks yield contributions to the scalar density that are similar to what would be 
obtained if the quarks are organized into nucleons. One part of the argument is based upon 
the well-known model-independent relation for the density dependence of the condensate 
[20] 

«p = ( 1 _J^ + ..A (3.10) 



(qq) V fl m l 

where is the pion-nucleon sigma term and p is the density of the matter. If we take 
U = 0.0942 GeV, = 0.138 GeV, p NM = (0.109 GeV) 3 and a N = 0.050 GeV, we find a 
reduction of the condensate in nuclear matter of 38%, which is consistent with relativistic 
models of nuclear matter [21, 22]. 

We now consider the corresponding relation for a quark gas of up and down quarks, 

+ ("I) 
(qq)o V fl m l J 

where p q is the density of quarks (p q = 3p) and a q is a "quark sigma term". We have shown 
in earlier work [23] that a q is in the range of 15-17 MeV, so that Eq. (3.10) and (3.11) imply 
that quite similar mean fields are generated by the quark gas and by nuclear matter. 

In Table I and in Fig. 2, we show the results obtained when Eqs. (3.8) and (3.9) are 
solved with G s = 9.00 GeV~ 2 , G D = -240 GeV" 5 , A 3 = 0.631 GeV, m° u = 0.0055 GeV and 
= 0.130 GeV. We note that the dependence of m u (p) on density is approximately linear 
for p/pnm < 2, with a 32% reduction in the value of m u (p) when pj ' pnm = 1- Another point 
to note is that m s (p) is density-dependent for finite values of Gp, since the (ss) condensate is 
modified by the coupling to the up and down quark condensates via the 't Hooft interaction. 
This coupling becomes less important as the up and down quark condensates are reduced 
at increasing density. [See Fig. 2.] 

We have also considered the solution for the SU(2) version of the above equations 

m u (p) =m° u -2G s (uu) p , (3.12) 

and have used the parameters specified in the Klevansky review article [10], Gs = 10.15 
GeV" 2 , m° u = 0.0055 GeV and A 3 = 0.631 GeV. The results for m u {p) are similar to that 
seen in Fig. 2, except that m u (0) = 0.336 GeV. [See Fig. 3.] In this case, m u (p) is reduced 
by about 32% when p = Pnm- 




FIG. 2: The solution of Eqs. (3.6) and (3.7) for the density-dependent constituent quark masses, 
m u (p) = m d (p) and m s (p) are shown. Here G s = 9.00 GeV~ 2 , G D = -240.0 GeV~ 5 , A 3 = 0.631 
GeV, ml = 0.0055 GeV and m° s = 0.130 GeV. 

IV. DENSITY AND TEMPERATURE DEPENDENCE OF THE CONFINING 
FIELD 

In part, our study has been stimulated by the results presented in Ref. [13] for the 
temperature-dependent potential, V(r), in the case dynamical quarks are present. We re- 
produce some of the results of that work in Fig. 4. There, the filled symbols represent the 
results for T/T c = 0.68, 0.80, 0.88 and 0.94 when dynamical quarks are present. This figure 
represents definite evidence of "string breaking" , since the force between the quarks appears 
to approach zero for r > 1 fm. This is not evidence for deconfinement, which is found for 
T = T c . Rather, it represents the creation of a second qq pair, so that one has two mesons 
after string breaking. Some clear evidence for string breaking at zero temperature and finite 
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L.3 
Hp 

( GeV 3 ) 


p/pnm 


m u (p) 
[GeV] 


m s (p) 
[GeV] 


0.00 


0.00 


0.358 


0.532 


0.007 


0.364 


0.318 


0.515 


0.010 


0.521 


0.300 


0.508 


0.0140 


0.729 


0.276 


0.498 


0.0192 


1.00 


0.242 


0.487 


0.025 


1.302 


0.200 


0.475 


0.030 


1.562 


0.162 


0.465 


0.035 


1.823 


0.121 


0.457 


0.040 


2.083 


0.0860 


0.452 


0.045 


2.343 


0.0618 


0.449 


0.050 


2.604 


0.0470 


0.448 


0.055 


2.864 


0.0378 


0.448 


0.060 


3.125 


0.0316 


0.448 


0.065 


3.385 


0.0272 


0.447 



TABLE I: Values of m u (p) and m s (p) obtained from the solution of Eqs. (3.6) and (3.7) are given 
for various values of the ratio p/pnm- (Here, kp = 0.0192 GeV 3 for nuclear matter, = 0.0055 
GeV, m° s = 0.130 GeV, A 3 = 0.631 GeV, G s = 9.00 GeV~ 2 , G D = -240.0 GeV" 5 .) 

density is reported in Ref. [24]. 

In order to study deconfinement in our generalized NJL model, we need to specify the 
interquark potential at finite density. We start with our model that was described in Section 
II. In that case we had V c (r) = /trexpf— fir). For the model we study in this work, we write 

V (r, p) = nrexp[— fx(p)r] (4.1) 



and put 



i-U- 



PC 



10 



350 




FIG. 3: The solution of Eq. (3.12) for m u {p) is shown. Here G s = 10.15 GeV~ 2 , m° = 0.0055 
GeV and A3 = 0.631 GeV. (See Table V of Ref. [10].) The dashed line is a linear approximation 
to the result which we use for p < 2p^M- (Nuclear matter density corresponds to k F = 0.0192 
GeV 3 . 

with p c = 2.25pnm and p = 0.010 GeV. With this modification our results for meson 
spectra in the vacuum are unchanged. Other forms than that given in Eqs. (4.1) and (4.2) 
may be used. However, in this work we limit our analysis to the model presented in these 
equations. The corresponding potentials for our model of Lorentz-vector confinement are 
shown in Fig. 5 for several values of p/pnm- 

We can see from Fig. 4 that, for T = 0.94T C , the use of dynamical quarks leads to 
an approximately constant value of V(r) = 1000 MeV for larger r. If we perform a Fierz 
rearrangement of the Lorentz-scalar potential to study pseudoscalar qq states, one introduces 
a factor of 1/4, making the value at large r to be about 250 MeV. (See Eq. (Bl) of Ref. [10].) 
However, rearranging the Lorentz-vector potential to study pseudoscalar qq states introduces 
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FIG. 4: A comparison of quenched (open symbols) and unquenched results (filled symbols) for the 
interquark potential at finite temperature [13]. The dotted line is the zero temperature quenched 
potential. Here, the symbols for T = 0.80T C [open triangle], T = 0.88T C [open circle], T = 0.80T C 
[open square], represent the quenched results. The results with dynamical fermions are given 
at T = 0.68T C [solid downward-pointing triangle], T = 0.80T C [solid upward-pointing triangle], 
T = 0.88T C [solid circle], and T = 0.94T C [solid square]. 

a factor of 1. Now, let us consider p/pnm = 0.94(pc/ pnm) — 2.11, and find the maximum 
of our Lorentz- vector potential at that density from the relation V max = k,/ p(p)e. Using our 
value for p(p) at p/pnm = 2.11, we obtain V max = 0.227 GeV. The value for the Lorentz- 
vector potential compares favorably with the value of V(r), for large r, quoted above. This 
result suggests that, if the dynamics of chiral symmetry restoration and deconfinement at 
finite temperature is somewhat analogous to the deconfinement process at finite density, our 
use of p/ Pnm = 2.25 may be a satisfactory choice. 
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1 2 3 4 5 



r (fm) 

FIG. 5: Values of V(r, p) are shown, where V(r,p) = nr exp[— p(p)r] and p(p) = po/[l — (p/pc) 2 ]- 
Here pc = 2.25pNM and po = 0.010 GeV. The values of p/pnm are 0.0 [solid line], 0.50 [dotted 
line], 1.0 [dash line], 1.50 [dash-dot line]. 1.75 [dash-dot-dot line], 2.0 [short-dash line], and 2.1 
[small dot line]. 

V. PION CONDENSATION AND THE CHOICE OF THE PARAMETERS OF 
THE INTERACTION 

It was suggested many years ago that the ground state of nuclear matter might have an 
unusual structure due to presence of pionlike excitations [25]. In finite nuclei such effects 
could imply anomalous behavior in states with J n = CT,l + ,2~..., etc. However, the 
nucleon-nucleon interaction is sufficiently repulsive in the relevant channel so that pion 
condensation does not take place at normal nuclear matter densities. That matter has been 
discussed in Ref. [26]. A constant g' parametrizes the strength of a nuclear force in the spin- 
isospin channel that represents short-range correlation effects and exchange effects. (See Eq. 
(5.11a) of Ref. [26].) The phenomenological value of g', obtained from the study of nuclear 
excitations, is sufficiently large so that pion condensation does not take place until about 
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three times nuclear matter density. (See Fig. 5.9 of Ref. [26].) 

In our work we will model the effects that prevent pion condensation by introducing a 
density-dependent interaction for the pionic states calculated in the NJL model. We write 

GM = G W (0)[1 - 0.087p/p NM ] , (5.1) 

where the second term in Eq. (5.1) represents medium effects that reduce the pion self- 
energy in matter. Here G^O) is the linear combination of Gs and Go given on page 269 of 
Ref. [12], 

G„ = G s + ^(ss). (5.2) 

Equation (5.1) represents our scheme for parametrizing the nuclear matter effects that pre- 
vent pion condensation. In our calculations of pionlike excitations we put G w (0) = 13.49 
GeV -2 , and used a constant values of G v = 11.46 GeV -2 . We may check that our choice of 
G w (0) is reasonable by using Eq. (5.2) with G s = 11.84 GeV" 2 and -180 GeV -5 < G D < 240 
GeV -5 . These values of Gs and Gd were obtained in our extensive study of the eta mesons 
[1]. Thus, if we take (ss) = -(0.258 GeV) 3 and G D = -190 GeV" 5 , we find G n (0) = 13.47 
GeV -2 . This analysis suggests that, once we fix our parameters in the study of the eta 
mesons, we can then infer the parameters needed for our study of the pion in vacuum. 

For this work, in our study of the kaon, we use G K (0) = 13.07 GeV -2 and G v = 11.46 
GeV -2 . Note that [12] 

G K (0) = G s + ( f(dd) o . (5.3) 

If we take G s = 11.84 GeV -2 , G D = -190 GeV -5 and (uu) = -(0.240 GeV) 3 , we find 
G K (0) = 13.15 GeV -2 , which is close to the value of G K (0) = 13.07 GeV -2 used in our 
calculations. In our work we have used 

G K {p) = G K (0)[1 - 0.087p/ PNM ] ■ (5.4) 

In the case of the kaon, about 40% of the assumed density dependence of Gk(p) may be 
attributed to the density dependence of (uu) p or (dd) p . We may consider the relation 

G K (p)=G s (p) + ^f(dd) p , (5.5) 
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and use a somewhat smaller reduction of Gs{p) for the kaon than that used for the pion in 
Eq. (5.5), since the reduction of (uu) p or (dd) p in matter effectively reduces the interaction 
strength. 

In the absence of a — f coupling we have G| 3 = G a() = G$ — (Gn/2){ss) [12]. If 
we again put G s = 11.84 GeV~ 2 , G D = -190 GeV" 5 , and (ss) = -(0.258 GeV) 3 , we have 
G a{) = 10.21 GeV -2 , which places the a (980) at 1.13 GeV. However, in the case of the scalar 
mesons there exist significant contributions to the interaction from processes that describe 
the scalar meson decay to two-meson channels. An extended discussion of these effects was 
given in an early work on scalar mesons [27]. In the case of the / (980) we presented a 
discussion of such terms as they affect the energy predicted for the /o(980) in Ref. [28]. 

In order to take into account these effects, which are not included in our RPA calculations, 
we increase the value of the ao coupling constant to G a() = 13.10 GeV~ 2 . That has the effect 
of moving the a (980) mass down to 980 MeV. 

We also introduce some density dependence of the interaction to avoid an "a condensate" , 
which would otherwise take place at p = 1.75pnm, if we use m u (p) = m^p) = 0.0055 + 
0.3585(1 — OAp/pnm)- Thus, we use G ao (p) = G ao (0)[l — O.OAdp/ pnm] when we allow for the 
rapid decrease in the value of m u (p) = m^p) given by the above expression. It is possible 
that the small reduction of G ao (p) in matter given above has it origin in a somewhat smaller 
attraction generated at the larger densities by the real part of the polarization operator that 
describes decay to the two- meson channels [27, 28]. We will provide further details of our 
treatment of the scalar mesons in Section VIII. 

VI. RANDOM PHASE APPROXIMATION FOR MESONIC EXCITATIONS 

In this work we report upon covariant random-phase-approximation (RPA) calculations 
of meson spectra in vacuum and in dense matter. Before writing the equations of our model, 
it is worth discussing some properties of RPA calculations made for many-body systems [29, 
30]. For example, such calculations have been performed to study excited states of nuclei. 
In the RPA one usually does not attempt to construct the wave function of the ground 
state. Rather, one considers amplitudes of particle-hole operators taken between the excited 
state and the ground state. The dominant amplitude usually involves the creation of a hole 
in the ground state and the creation of a particle in what are predominantly unoccupied 
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states. Smaller amplitudes are found if one destroys a hole in the ground state and destroys a 
particle in the predominantly unoccupied states. These smaller amplitudes are only nonzero, 
if one allows for correlations in the ground state. 

Such RPA calculations are particularly important for states that are collective with re- 
spect to matrix elements of electromagnetic transition operators, for example. In hadron 
physics the most "collective state" is the 7r(138). In this case the "large" and "small" com- 
ponents of the wave function, in the sense of the RPA, are comparable in magnitude and 
approach equality in magnitude as one approaches the chiral limit, when m n — > 0. 

Another important feature of RPA calculations is that they may be considered as an 
investigation of the properties of small oscillations about the ground state. Thus, if one 
obtains an imaginary energy value for the ground state, one infers that the ground state is 
unstable. A new ground state must be constructed that will yield real eigenvalues. (Note 
that imaginary eigenvalues may be obtained, since the RPA Hamiltonian is not Hermitian.) 

There is a strong analogy that can be made between the particle-hole RPA calculations 
described above and the calculation of mesonic excitations. For example, a "hole" in the 
ground state (the vacuum) is an antiquark, while the particle state is the quark. If we 
perform relativistic RPA calculations for the pion and its radial excitations, an imaginary 
energy calculated for the pion is a signal of pion condensation. 

Random-phase-approximation equations may be derived using the vertex equation of Fig. 
lb. The RPA equations for the study of the pion, kaon, and eta mesons were derived in Ref. 
[1]. In the case of the pion and kaon we include pseudoscalar — axial- vector coupling. The 
most complex case is that of the eta mesons which, in addition to pseudoscalar — axial- vector 
coupling, involves singlet-octet coupling in the flavor sector. 

In this work we only record the equations in the simplest example, that of RPA calcu- 
lations for the a mesons [31]. In this case the large component is denoted as 4> + (k), while 
the small component is (p~(k). These functions are found to satisfy coupled equations for 
mesons in vacuum: 



2E u (k)<f ) + (k) + / dk' [H c (k,k') + H NJL (k,k')](f> + (k') 



(6.1) 
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(6.2) 



where E u (k) = [k 2 + m 2 u ] 1/2 , 



H c (k,k') 



1 [2V c (k,k')k 2 k' 2 + kk'V 1 c (k,k')] 
(2tt) 2 E u (k)E u (k>) 



(6.3) 



and 



8N C k 2 k' 2 G ao e- k2 l 2a2 e- k ' 2 ' 2a2 
(2tt) 2 E u {k)E u {k') 



(6.4) 



In Eq. (6.3) we have introduced 



v l c {Kk>) = 1 -]_ 



-i 



dx Pi(x)V c \k - %') . 



(6.5) 



Here, x = cos9 and Pj(x) is a Legendre function. The terms exp[— k 2 /2a 2 ] and exp[— k' 2 /2a 2 } 
are regulators with a = 0.605 GeV. 

In order to solve these equations in the presence of matter, we replace m u , G ao and 
p by m u (p), G ao (p) and p(p). (Recall that p(p) = po/[l — (p/Pc) 2 ]-) 111 our calculation 
for the a states we have taken m u (p) = m" + 0.3585 GeV [1 — 0Ap/p NM ] and G ao (p) = 
G ao (0)[l — 0.045p/pjvm], with m° = 0.0055 GeV. As an alternative, the mass values for 
m u{p) — nT>d(p) m ay be taken from Table I. 

VII. RESULTS OF NUMERICAL CALCULATIONS: PSEUDOSCALAR MESONS 

The choice of the parameters in the case of the pion and its radial excitations was discussed 
in Section V. We use G n (p) = G 7r (0)[l - 0.087p/pAr Af ] and m u (p) = m d (p) = 0.0055 + 
0.3585[1 - 0Ap/p NM ] with G ff (0) = 13.49 GeV -2 and G v = 11.46 GeV" 2 . Also, p(p) = 
Po/[l — (p/pc) 2 ] with p = 0.010 GeV and p c = 2.25p NM - 

The results of our calculations are shown in Fig. 6. At p = 0, the first radial excitation 
of the pion is found at 1.319 GeV. The large number of states above 1.3 GeV have wave 
functions that are dominated by either the 75 or 7075 vertex. The pion wave function has 
mainly a 75 vertex structure, with a small admixture of the 7075 vertex. (The axial-vector 
part of the wave function makes a significant contribution in the calculation of the pion 
decay constant, f n .) 
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FIG. 6: The mass values for the pion and its radial excitations are presented as a function of 
the density of matter. Here, G n {p) = G W (0)[1 — 0.087p/ pnm] and m u (p) = m^{p) = + 
0.3585 GeV[l - 0Ap/p NM ], with m° u = 0.0055 GeV. We use G n (0) = 13.49 GeV~ 2 , G v = 11.46 
GeV~ 2 and p = p /[l - (p/p c ) 2 }, with = 0.010 GeV and p c = 2.25p NM - 

It may be seen from the figure, that with the reduction of the value of the constituent 
mass and of the confining field with increasing values of p/pnm, the radial excitations that 
appear as bound states become fewer in number. Beyond p/pnm = 1-50 only the nodeless 
pion wave function is bound and that state is no longer supported beyond p/pnm — 1-80. 
That represents the beginning of the deconfined phase in the case of the pion for the model 
introduced in this work. 

Somewhat similar behavior is found for the kaon and its radial excitations, as may be 
seen in Fig. 7. Here we have used the mass values given in Table I and G K (p) = G^(0)[1 — 
0.0S7p/p NM ] with G K (0) = 13.07 GeV~ 2 and G v = 11.46 GeV -2 . Again we see only a 
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FIG. 7: Mass values of the K mesons are shown as a function of the density of matter. Here 
we use G K (0) = 13.07 GeV~ 2 , G K (p) = G K (0)[1 - 0.087p/p NM ], G v = 11.46 GeV" 2 and p = 
po/[l — {p/pc) 2 ], with po = 0.010 GeV and pc = 2.25pNM- The mass values given in Table I are 
used. 

small increase of the mass of the nodeless state, the pseudo Goldstone boson, as p/pnm 
is increased. We again find deconfinement for p/pnm > 1-8. The density dependence of 
Gk(p) is taken to be the same as in the case of the pion. However, in this case, we have 
noted previously that about 40% of the reduction of Gk(p) with increasing density may be 
ascribed to the density dependence of the up and down quark condensates, (uu) p and (dd) p . 
The calculation of the density dependence of the coupling constants in our model is a major 
undertaking and is beyond the scope of this work. 
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VIII. RESULTS OF NUMERICAL CALCULATIONS: SCALAR MESONS 



We have recently discussed the properties of the / (980), giving particular attention to 
the role of the polarization diagrams that describe the decay of the f mesons to the nn 
or KK channels [28]. (See Fig. 2 of Ref. [28].) However, when we diagonalize the RPA 
Hamiltonian we do not take those terms into account. Calculations of such effects are more 
easily made if we construct a quark-antiquark T matrix. For a single channel example we 
may write 

^--r^wr (81) 

where G is the appropriate coupling constant for that channel and J(p 2 ) is the corresponding 
vacuum polarization operator. In our model J{p 2 ) is calculated with the confining vertex 
function that appears in Fig. la as a crosshatched region. (See Fig. 1 of Ref. [28].) The 
resulting J(p 2 ) is a real function, which is singular at the values of p 2 for which there is a 
bound state in the confining field. If we include polarization diagrams that describe coupling 
to two-meson decay channels, Eq. (8.1) is modified to read 

tijP ' = ~ 1 - G[J{p 2 ) + ReK(p 2 ) + ilmK{p 2 )} ' (8 ' 2) 
The calculation of J(p 2 ) and K(p 2 ) has been extensively discussed in our earlier work. In 
the case of the scalar mesons, inclusion of ReK(p 2 ) can move the mass of the lowest-energy 
state down by about 70-100 MeV [27, 28]. 

In the case of the ao(980), the use of Gs and Gd determined in our study of the eta 
mesons places the ao(980) at 1.13 GeV. In the present work we have increased the coupling 
constant from G ao = 10.21 GeV -2 to G ao = 13.10 GeV -2 to move the lowest a state 
down to 980 MeV. That creates a problem of "a condensation" which we avoid by taking 
G ao (p) — G ao (0)[l — 0.045p/ Pnm}- One may speculate that the effects that increase the 
effective coupling strength from G ao = 10.21 GeV~ 2 to G ao = 13.10 GeV~ 2 have some 
density dependence that reduces the induced attraction at the higher densities. 

In Fig. 8 we show our results for the a$ mesons. There we see deconfinement at about 
p = 2.0p7VM which is a slightly larger value of the density than that found for the other 
mesons studied in this work. However, the behavior of the lowest ao state with increasing 
density is made somewhat uncertain because of our lack of knowledge of the appropriate 
form for G ao (p). 
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FIG. 8: Mass values for the ao mesons are given as a function of the matter density. Here, 
we have used G ao (0) = 13.10 GeV" 2 and G ao (p) = G ao (0)[l - 0M5p/p NM }- We have used 
m u = m o + 0.3585 GeV[l - 0Ap/p NM ] with m° = 0.0055 GeV. The dotted line results, if we put 
G ao (p) = G ao (0)[l — 0.087 p/pnm] and use the mass values of Table I. The dotted curve is similar 
to the curve for the ao mass given in Ref. [31]. The curves representing the masses of the radial 
excitations are changed very little when we use the second form for G ao (p) given above. 

For our study of the f mesons we work in a singlet-octet representation and use the 
coupling constants = 14.25 GeV -2 , Gf 8 = 10.65 GeV~ 2 and = Gf = 0.4953 
GeV" 2 . This choice yields 980 MeV for the mass of the / (980). The fact that > G§ 8 
is a feature of the 't Hooft interaction and leads to the /o(980) being mainly a singlet state 
[28]. (For the r?(547) the behavior of the 't Hooft interaction is such that Gf 8 > Gg, [12, 28] 
and, therefore, the ??(547) is predominantly a flavor octet meson [1].) 

In our study of the f mesons at finite density we use the mass values of Table I and do 
not introduce any density dependence for Gq , Gf 8 and G^g. The results of our calculation 
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3.5 




FIG. 9: The figure shows the mass values of the /o mesons as a function of density. The mass 
values for the quarks are taken from Table I. In a singlet-octet representation, we have used the 
constants = 14.25 GeV" 2 , = 10.65 GeV" 2 and G§ 8 = 0.4953 GeV~ 2 . Deconfinement 
takes place somewhat above p = l.8pNM- Here p = po/[l — (p/pc) 2 } with po = 0.010 GeV and 
pc = 2.25p NM - 

are shown in Fig. 9. Since the /o(980) has a significant ss component, the mass value only 
decreases slowly, with a value of 700 MeV for the lowest f state at p/pnm = 1-82, where 
deconfinement sets in. 

In Ref. [28] we provide a discussion of the T matrix for the singlet-octet channels. There 
the role of Kq (p 2 ), Kq 8 (p 2 ) and K 88 (p 2 ) in lowering the energy predicted for the / (980) is 
discussed in some detail. 

Our results for the energy levels of the Kq mesons are given in Fig. 10. In this case we 
use a constant value for G K * = 10.25 GeV -2 . The results are hardly modified if we allow 
for a small density dependence of Gk*- Since the Kq mesons contain a strange quark, the 
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FIG. 10: The figure shows the mass values obtained for the Kq mesons as a function of density. Here 
we use a constant Gk* = 10.25 GeV -2 . Deconfinement takes place somewhat above p = 1.8pjvAf- 
The quark mass values were taken from Table I. 

density dependence of their energies is not as marked as that of the ao mesons which only 
contain up and down quarks in our model. In that regard, the behavior of the Kq mesons 
is more like that of the fo mesons, which have some strange quark content. Again we see 
deconfinement for p > 1.8pnm- 

IX. DISCUSSION 

We originally chose pc = 2.25p NMl since the curve in Fig. 2 that shows the values of 
m u {p) seemed to change its behavior at about k% = 0.045 GeV 3 , which corresponds to 
p ~ 2.3p NM - We can attempt to see if that is a reasonable choice by noting that "string 
breaking" should occur when the energy of the extended string is equal to the energy of 



23 



the lowest two-meson state that can be formed when the string breaks. Therefore, we may 
write V m ax = mi + m 2 , where mi + m 2 are the masses of the mesons in the final state. 
We then use V max = K,/p(p)e to find a value p{p) and obtain p/pc from the expression 
p(p) = po/[l — (p/pc) 2 ]- We then put pc = 2.25p NM and calculate the value of p/pnm 
where we might expect string breaking. We consider the final states 7T7t, nK, ni] and KK. 
The corresponding values of p/ Pnm are 2.09, 1.86, 1.83, and 1.61 for p c = 2.25pnm- Note 
that the fT(495) and ^(1430) mesons can break up into the irK system, while the a (980) 
is strongly coupled to the irr/ channel. The / (980) is coupled both to the irir and KK 
channels. On the whole, the values of p/pnm calculated above are generally consistent with 
the value of that quantity that leads to deconfinement in our model. That result tends 
to suggest that, for light mesons, the density that leads to string breaking may be similar 
to the density for deconfinement. (In general, however, these processes are distinct and 
further studies would be needed to see if string breaking and deconfinement are related at 
finite density.) We may suggest that, if the initial meson is of the same type as the mesons 
that appear upon string breaking, it becomes reasonable to suggest that the instability of 
the initial mesons is also felt by the final state mesons, giving rise to the relation of string 
breaking and deconfinement suggested above for light mesons. 

A comprehensive discussion of meson properties at finite temperature and density has 
been presented by Lutz, Klimt and Weise [32]. Since those authors did not include a model 
of confinement, they were able to calculate values of the meson masses for large values of the 
density. Their Fig. 8 shows the calculated masses of the nodeless pion, f and a mesons 
for < p/ Pnm < 3.5. They also give the result for an f' Q excitation. (The f and f exhibit 
singlet-octet mixing.) Compared to our results, their value of the fo mass falls more rapidly 
than ours, becoming degenerate with the pion mass at about p/pnm = 3. On the other 
hand, the mass of the a in their work is about 600 MeV at pj ' Pnm = 2. They are able 
to derive systematic low-density expansions for various quantities which provide important 
insight into the results obtained in numerical studies. They also show that effects due to 
finite quasiparticle size are important in stabilizing the density and temperature dependence 
of the pion mass. The main deficiency of their work is the absence of a model of confinement. 
Therefore, we believe our work provides a natural extension of the work reported in Ref. 
[32]. 

It is worth noting that deconfinement takes place in our model at about p = 1.8pnm, 
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while the confining potential goes to zero at p = p c = 2.2hpNM- That suggests that the 
specific form we have chosen for the density dependence, p(p) = fJ>o/[l — (p/pc) 2 }, is not 
particularly important. What is more important is the behavior of our confining potential, 
V (r, p), shown in Fig. 5. There, we see that the potential still has a substantial magnitude 
at p = l.75p NM and p = 2A0p NM . 

Since the analysis of Ref. [32] is made in the absence of a model of confinement, many 
analytic results can be obtained for the behavior of various quantities when small changes 
in density and temperature are considered. Indeed, the work of that reference provides 
some support for our treatment of the pion and kaon. It is shown that the Goldstone boson 
remains at zero mass in the chiral limit as long as the system remains in the Goldstone- 
Nambu mode of symmetry breaking. For finite current quark masses, we quote the result 
given in Eq. (5.6) of Ref. [32] for T = 0, 

dml , d(uu) 

f= {l-2m 2 u (r 2 )) JL> . (9.1) 



mi / ?n ^ 



Here, rs is the quasiparticle radius. That quantity is defined in terms of the form factor 
F s (p — p') that appears in the matrix element of the u-quark scalar density 

(u(p')\uu(0)\u(p)) =Fs(p-p')u(p')u(p). (9.2) 

In Eq. (9.2) u{p) denotes the Dirac spinor of a constituent u quark with four-momentum p. 
The scalar mean-squared radius is then 



(r 2 s) = ^nF s ( q 2 : 



(9.3) 

q 2 =0 



(See Eq. (A. 7) of Ref. [32] for an explicit expression for (r|) in terms of the parameters of 
the NJL model.) With the well-known relation [20] 

d(uu) o^p 



(uu) mlfl ' 



(9.4) 



Eq. (9.1) becomes 



dml = -(l-2ml(rl))^f. (9.5) 

J 7T 



If one ignores the quasiparticle size, one has drm\ = —{a^p/f 2 ) [33, 34], which implies pion 
condensation at a critical density p cr u = flrm\/a^ = (0.148 GeV) 3 , which is about 2.5 Pnm ■ 
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The second term in Eq. (9.5) works against condensation. With m u = 0.364 GeV and 
rs = 0.40 fm [32] one finds that 6m% increases slowly with increasing density, as born out 
by the calculations reported in Ref. [32]. Our choice of G w (p) = G 7r (0)[l — 0.087p/ Pnm] 
reproduces the almost constant value of m^. We see that the density-dependent term in 
G v (p) plays a similar role in our model as that played by the second term in Eq. (9.5). 

We have some confidence in our treatment of the pion and kaon at finite density. We 
recall that we were able to find satisfactory values of G„(0) and Gk(0) using the values of Gs 
and Go obtained in our study of the eta mesons [1]. Therefore, our work provides a unified 
approach for the nonet of pseudoscalar mesons in the presence of a model of confinement. 

Since confinement is important for the ao(980) and /o(980) mesons, it is uncertain whether 
the results of Ref. [32] for the properties of these mesons can be trusted. These mesons are 
in the continuum of the NJL model without confinement and various assumptions need to 
be made as to how the formalism is to be applied. For a small increase in density, the mass 
of the a in our model and in Ref. [32] are similar. For the larger values of density, the use 
of G ao (p) — G ao (0)[l — 0.0A5p/pnm] in model leads to a rather small mass for the a for 
p ~ 2p NM . [See Fig. 8.] 

Our treatment of the ao mesons is less satisfactory than that of ir and K mesons, since 
coupled channel effects are important in the case of the scalar mesons. Using the values of 
Gs and Go obtained in our study of the eta mesons [1], we found the lowest a state at 1.13 
GeV. To place the a at 980 MeV, we increased the value of G ao (0). That increase led to 
the possibility of an ao condensation, which was removed by reducing the coupling constant 
with increasing density. [See Fig. 8.] However, it might be preferable to accept the value of 
1.13 GeV for the mass of the ao and, therefore, avoid the problem of ao condensation. Our 
difficulty in this case arises since we do not know the density dependence of the processes 
that move the a mass from our predicted value to the experimental value of 980 MeV. 

In our model we see some relation between the partial restoration of chiral symmetry 
and deconfinement. With reference to Fig. 2, we see that the up (or down) quark mass 
drops in a roughly linear manner with increasing density up to about 2 or 2.5 times nuclear 
matter density. With the reduction of the magnitude of the confining field, as seen in Fig. 
5, the combined effect of the smaller confining field and reduced quark mass values leads to 
deconfinement at about 1.8 Pnm ■ For the fo, K and Kq mesons, the reduction of the mass 
of the quarks is less important, since these mesons have either one strange quark (K, Kq) 
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or an ss component (/o). However, deconfinement still takes place at about p = l.Sp^M for 
the mesons. 

In future work we will study the dependence of the deconfinement process on both tem- 
perature and density. In addition, it would be desirable to have some understanding of the 
mechanism by which the increased matter density modifies the confining interaction. 
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